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Abstract 

In hadronic collisions, the mini-jet cross section is formally divergent in the 
limit p± ^ 0. We argue that this divergence is tamed by some effective 
colour correlation length scale of the hadron. A toy model of the hadronic 
structure is introduced, that allows an estimate of the screening effects, and 
especially their energy dependence. 



The perturbative parton-parton interaction cross section in a hadronic collision is 
divergent roughly like 

da _ 



oc 



dpi p1 

for p_L — i> 0. This is to be convoluted with parton distributions fi{x,Q'^ ~ p^) to give a 
hadronic interaction cross section. The dominant contributions come from scattering by 
t-channel gluon exchange. With an artificial lower cut-off scale p±min, 

a-.nt{p±m:n) = E / Q') J dx, /.(x^, Q') dpi "^^'^^'^^^ = ^^^^"^ . (2) 
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The jet/mini-jet cross section is twice this, since each scattering gives two jets in the 
lowest-order approximation considered here. Studying e.g. Tevatron collider energies, 
y/s ^ 2 TeV, o"int(p±min) exceeds the total pp cross section cxtot for p±_mm ~ 3 GeV. Since 
3 GeV still is well above the Aqcd scale of ~ 0.2 GeV, there is no obvious reason why 
perturbation theory should have broken down, and one would seem to be in trouble. The 
resolution of this paradox probably comes in several steps, and at least in two. 

First of all, the mini-jet cross section above is inclusive. Thus, if an event contains 
two parton-parton interactions it counts twice in o"int(p±min) but only once in atot- 
Thereby crint(p±min) > ^tot becomes allowed. Multiple parton-parton interactions is the 
concept that, based on the composite nature of hadrons, indeed several parton pairs 
may scatter in a typical hadron-hadron collision 0. Over the years, evidence for this 
mechanism has accumulated 0, such as the direct observation by CDF ||^. The events 
studied experimentally, with two parton pairs at reasonably large p±, only form the tip of 
the iceberg, however. One may expect that most interactions are at lower where they 
do not produce visible jets, but only contribute to the underlying event structure. As 
such, they are then believed to be at the origin of a number of key features, like the broad 
multiplicity distributions, the significant forward-backward multiplicity correlations, and 
the pedestal effect under jets 0. 

While the effects of multiple interactions on event properties are smaller for lower- 
p± scatterings, in the models we studied so far the interaction rate increases faster with 
decreasing p± than the effect per scattering, so there is no stability in the limit p±min — ^ 0. 

The second necessary aspect is likely that a regularization of the jet cross section 
should occur at small p± from the fact that the incoming hadrons are colour singlets — 
unlike the coloured partons assumed in the divergent perturbative calculations — and 
that therefore the colour charges should screen each other in the p± ^ limit. Thus 
P±min could take on a physics meaning, roughly, as the inverse of some colour screening 
length in the hadron. Of course, one would not expect a sharp cut-off of the mini-jet 
cross section at p_Lmin, but rather a smooth dampening. Nevertheless the p_Lmm parameter 
provides a useful first approximation. Fits to data typically give p_Lmin ~ 2 GeV. 

One key issue is the energy dependence of p±min] this may be relevant e.g. for compar- 
isons of jet rates at different Tevatron energies, and even more for any extrapolation to 
LHC energies. The question actually is more pressing now than at the time of the study 
in 0, since nowadays parton distributions are known to be rising more steeply at small x 
than the fiat xf{x) behaviour normally assumed for small before HERA. This trans- 
lates into a more dramatic energy dependence of the multiple-interactions rate for a fixed 
P±min, and thence to a charged multiplicity rising faster than data p. Based on such con- 
siderations, the p_Lmin in the Pythia program was made explicitly energy-dependent 
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some time ago: 

P±min = (1.9 GeV) (p^y , (3) 

with e = 0.08. This value is picked to agree with the s"^ behaviour assumed in the 
parameterization of the total cross section in hadron-hadron collisions ^ that, via reggeon 
phenomenology, should relate to the behaviour of parton distributions at small x and Q^. 
A study of the energy dependence of dNcharged/ (iv\v=o ^^^^ confirms that an e in this range 
gives sensible agreement with data [§]. 

In the following, we will develop an extremely simple model for the proton structure, 
where the amount of screening may be studied 0. Our basic approach is to assign 
an exclusive parton configuration to each proton. This configuration is — for a given 
evolution scale Q — defined in terms of the longitudinal momentum fraction x values, the 
transverse coordinate positions and colour charges of 'all' the partons present in a proton. 
The traditional inclusive parton-distribution picture is obtained as the average over many 
such proton 'snapshots'. A simple plane wave represents the exchange of a gluon between 
two colliding hadrons. Longitudinal coordinates of the partons are not specified, so the 
exchanged gluon is only characterized by its transverse momentum. The exclusive proton 
picture then allows a comparison of a coherent versus an incoherent sum of colour charges 
coupling to the plane wave. It is this ratio, coherent over incoherent, that will provide our 
measure of screening between partons in the proton. In particular, the long wavelength 
limit then explicitly corresponds to a vanishing of interactions between two colourless 
hadrons. As the collision energy is increased, more partons at smaller x values become 
accessible and the screening effects increase in importance. 

We certainly know this ansatz to be wrong in many of its details, but still hope that 
it will catch enough of the spirit to provide some insight. It is clear that a full-fledged 
description — e.g. with lattice QCD calculations — is way beyond the current capability, 
so toy models is all that could be offered today. The approach we present here should be 
viewed as only one possible aspect or formulation of the required dampening effect. The 
vanishing net colour charge of hadrons has been used in models since long |jlO[, as has the 



spatial size of the colour singlet wave function W^. Our model also has some similarities 



with the dense-packing view of the proton |12|, although without the possibility of parton 



recombinations. Among alternatives, one could mention the possibility of a nontrivial 



structure of the QCD vacuum |[T^, the solution of a Dyson-Schwinger equation for the 
gluon propagator [jl4|, the introduction of an effective gluon mass from lattice QCD results 
T5| , and the possibility of several hard scatterings within a single parton chain [|16], e.g. 



formulated in terms of non-integrated structure functions and the Linked Dipole Model 
description of non-p^-ordered chains |]17 . 



In our model, the evolution of the parton distributions is given by the conventional 
DGLAP equations As usual, the starting configuration requires nonperturbative 



input. A valence- like ansatz at a small Qo scale [|T9| limits the number of required free 
parameters. We have therefore used the GRV distributions as a reference [^, where 
the three valence quarks together with two 'valence gluons' almost completely define the 
proton at the small scale Qo = 0.48 GeV. In our simulation, these five partons are first 
selected according to a single-particle distribution 

/(x) ocx"(l-a;)^ , (4) 

with tuned parameters a = —0.4 and (3 = 1.2. Thereafter all five Xi are normalized so 
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Figure 1: Parton momentum distribution, summed over quarks and gluons. Our sim- 
ulation is compared with the GRV parameterization pO[. (a) At the small scale 
Qq = 0.48 GeV. (b) Evolved to = 10 GeV^, with a cut so that branchings do not 
produce partons below Xmin = 10~^. 



that their sum corresponds to unity: 

(3;i)norm = :^^5 ■ (5) 

This means that the original distribution is pushed towards the middle, and it is especially 
difficult to tune the low- a; tail, at least with the simple f{x) ansatz above. A reduction 
of the Qq scale from 0.48 to 0.44 GeV helps improve the small- a; behaviour. Fig. |l|a. 

From there on, the evolution follows the standard leading-order DGLAP formalism, 
expressed in terms of Monte Carlo-generated branchings of partons, q qg, g — gg and 
g — i> qq. Unfettered, such an approach gives an infinity of gluons in the x ^ limit. 
For the exclusive approach we have in mind, this is not so physical. In the study of 
scatterings above some p±ram scale, only partons above some related Xmin scale can at all 
interact. Partons below Xmin then have a dubious existence. To simplify the picture, we 
assume they are not resolved, i.e. that only branchings producing both daughters above 
Xmin are allowed. This reduces the amount of evolution, but not as dramatically as might 
be guessed at first glance. In Fig. we see that, in a typical case, the parton density 
with an x^i^ constraint is at most 20% above the no-constraints normal evolution in the 
region x > x^i^, while the former obviously vanishes for x < a;min- Disagreements are 
more visible — but also less interesting for us — in individual parton species, especially 
for sea quarks, which are not part of our simple ansatz at Qo but are in the GRV one. 

The kinematics of a hard scattering requires that xiX2 = s/s > 4p\/s. For collisions 
at central rapidities this suggests Xmin — 2p_L/i?cM- Away from the center, Xmin would be 
smaller on one side, down to Ap^/E^-^, and correspondingly larger on the other, which 
to some extent should cancel out. Therefore we will have as standard scenario that Xmin 
should increase proportionately to the p± considered, but will include as alternative a fixed 
Xmin- Furthermore, we will make the association that Eqm oc 1/xmin, for Xmm evaluated 
at some fixed reference pj_ like 1 GeV, but allow as alternative that the relation could be 
more like Eqm oc l/^/x^. 

Next we consider the assignment of transverse position coordinates. The original 
five partons are here selected according to a Gaussian shape, with a radius of 0.7 fm 
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in each of the two transverse dimensions. This gives a minimal correlation between the 
position coordinates of these five partons, and thus probably errs on the side of simplicity. 
Attempts with somewhat more correlated forms gave fairly similar results, however, so 
the Gaussian ansatz is not so critical. 

The spatial extent of a fluctuation like q — > qg — > q, g — > gg — > g and g — > qq — > 
g is of order IjQ according to the uncertainty relation. Therefore, when a DGLAP 
branching occurs at a scale Q, the two daughters are assumed to have time to fluctuate 
a distance of this order away from their production coordinates before they branch in 
their turn, or are probed by our plane gluon wave. The daughter partons are assumed 
to move out in opposite directions from the production vertex, with a uniform azimuthal 
distribution. The two distances are picked independently of each other, uniformly between 
and l/Q, alternatively between and 2/(5. We will see that this choice is relevant, 
whereas results are essentially equivalent e.g. for exponentially dampened distributions 
with the same mean. For a parton which takes the momentum fraction z in a branching, 
an additional factor \/\ — z is included, giving an approximate pj^/Q dampening factor 
for the separation in transverse rather than longitudinal direction. Here the relation 
P± = Vl — zQ follows in light-cone kinematics assuming the mother and recoiling partons 
to be massless. It is of some relevance that a very soft gluon emission does not affect the 
position of the hard parton unduly, since else multiple soft gluons could lead to a too 
rapid increase of the proton radius with Q. For the final result here, however, the yjl — z 
factor turns out not to be crucial. 

Finally we come to the colour space picture. Here we have picked a simple planar 
representation, with the three primary colours placed in a triangle around the origin. 
Thus the full phase information is lost; i.e. it is possible to ensure that the proton state 
is colour neutral but not that it is in a singlet. This certainly is a major simplification of 
the real world, and one which it is difficult to estimate the impact of. A gluon emission 
corresponds to a quark changing colour, e.g. q(r) — > q(&) + g{rh). The emitted gluon 
thus carries one colour and a different anticolour, i.e. the two colour diagonal gluon states 
are not populated; another simplification. Further, the ratio of gluon to quark charge is 
v^S versus \J Nc/Cp = 3/2 in QCD. Also in a gluon branching to two, a new colour is 
picked, e.g. g{rb) —>■ g{rg) + g{gb) , in order to avoid diagonal gluons. The gluon branching 
to quarks is unambiguous, g{rb) — > q(r) +q(b). For colour assignments in the starting 
configuration, the three valence quarks are first picked to be r + + 6, and then the two 
gluons arc handled as if emitted from two of the quarks, picked at random. 

An effort is thus made to retain overall vanishing colour for all partons in a proton that 
are seen by an exchanged gluon. In case of multiple interactions, it could be argued that 
only the 'first' exchange would be between singlet hadrons, while 'subsequent' ones would 
be between already coloured hadrons. Clearly, the imposition of such a time ordering 
can be questioned but, even if taken at face value, colour screening could still provide a 
significant dampening of the naive perturbative answer. 

To summarize, an explicit parton configuration of the proton is constructed for any 
value of the evolution parameter Q above some low Qq starting scale. Each parton k 
is characterized by its momentum fraction Xk, its transverse position coordinate and 
its two-dimensional colour charge qk- If such a proton is probed by a gluon plane wave 
function of transverse momentum p^, we may define a ratio 
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Figure 2: (a-c) The amount of screening, A in eq. as a function of the probing gluon 
(transverse) momentum P = p±. (a) For a proton at fixed evolution scale = 10 GeV^, 
and fixed Xmin = 10~^. (b) For a proton where the proton evolution scale Q = p± for 
P± > 0.44 GeV, alternatively 1 GeV, again with Xmin = 10"'^. (c) Full: variable Xmin = 
10~^p_L and fixed Q = 2 GeV. Dashed: variable = lO^'^px and Q = max(p_L, 1 GeV). 
Dotted: fix Xmin = 10~^ and variable Q = max(p_L, 1 GeV). (d) The Xmin dependence of 
the pj_o effective cut-off. See text for a explanation of the three scenarios. 



between a coherent and an incoherent sum of the colour charges. Note that we assume 
the exchanged gluon to be a superposition of all possible colour-anticolour pairs, so that 
we do not need to single out a special direction in colour space. 

The screening effect is most easily observed if one fixes the Q of the probed proton 
and then varies the p± of the probing gluon. Fig. |^a. The suppression factor A then 
has to vanish quadratically a.s p± — > 0, but in the interesting p± range the dominant 
aspect is instead the slower rise to the asymptotic limit A = 1. The choice of a fixed 
Q scale is representative for the way parton distributions de facto are used with Qq as 
scale for Q < Qq. For popular distributions like the current CTEQ and MRST ones pi 



where Ql > 1 GeV , the inclusion of a suppression factor then significantly modifies the 
picture for Q < Qq. Actually Qq has tended to come down with time; some years ago 
Qq > 4 GeV^ was the norm. Also remember that the cross section suppression is A^, i.e. 
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with one factor of A for each hadron beam. 

A conventional choice, but not a unique one, is to pick the Q scale of parton distribu- 
tions to agree with the p± scale of the hard interaction. Then the approach A — > 1 is much 
slower at large p±, since the branchings that occurred at scales Q not so much smaller than 
p± are not fully resolved. Fig. We remind that the proton seen by a gluon probe is a 
much more 'busy' place than the one seen by a photon in Deeply Inelastic Scattering pro- 
cesses, where the gluon content is not probed directly. Thus the destructive interference 
we introduce here does not invalidate the conventional structure function measurements 
and partonic interpretations. In hadronic physics, higher order perturbative corrections 
could partly contain, and partly mask, the physics of an A somewhat below unity. 

More interesting is the behaviour at small p±. We here note a spike in A at the scale 
P± = Qo = 0.44 GeV, Fig. At this scale, the proton consists of the original five partons, 
and these are normally well separated in space. Above it, the partons emitted in the 
evolution process tend to have a more clustered spatial distribution and therefore screen 
more. Below it, there is no further evolution of the proton but only a straightforward 
destructive interference. Thus the spike is an artifact of how we discontinuously in Q go 
from a fixed proton picture to a very rapid evolution; remember that the evolution rate 
behaves like Q!s(Q^) dQ^/Q^ and thus is as largest just at the lower cut-off. In order to 
provide a smoother and physically maybe more sensible picture, we therefore do not probe 
the proton at scales below 1 GeV, i.e. put Q = ma.x{p±, 1 GeV). The original Qq scale 
and evolution below 1 GeV still lives on in terms of a sensible x/r /colour configuration 
at 1 GeV. 

Fig. 0b is based on a fixed Xmin scale. However, an Xmin scaling oc p± behaves very 
similarly. Fig. if the two are matched for p± = 1 GeV. If instead Q is fixed, the shape 
is more changed. 

We now turn to the prime objective of this paper, namely to study the energy depen- 
dence of p±mm- Actually, the sharp cut-off 0{p± — p_Lmin) is not very physical. A more 
likely behaviour is a dampening factor something like p\/ {p\ +p'±o), where p±Q takes over 
the role of free parameter. With one such factor per incoming proton, the l/p^ singularity 
of the perturbative parton-parton cross section is regularized. 

Whereas the proposed dampening form does not completely reproduce the curves 
shown above, qualitatively there is agreement. Instead of doing a fit to the curve shape, 
p_Lo is extracted as the p± value for which A = 1/2. This would also have been a sensible 
definition of a p±min, backed up by experience with the impact of the two regularization 
procedures on fits to data, where pj_o ~ p±min is obtained . 

The resulting p±o values are plotted as a function of 1/xmin in Fig. |^. Three scenarios 
are compared: 

1. The proton probed at a fixed Q = 2 GeV scale, travel distance of a parton is 
distributed between and 1/Q of the branching where it is produced. 

2. Also proton probed at a fixed Q = 2 GeV scale, but travel distance twice as large. 

3. The proton is probed at a running Q = pj_, travel distance between and 2/Q as 
in 2. 

In all these three cases, representative for several more studied by us, a clear energy- 
dependence of pj_Q is observed. If one tries to make a fit to an exponential form pj_q oc 
l/^miii) typically one obtains a 5 in the range 0.05 to 0.08. However, two comments are 
in order. One is that the curves tend to rise less steeply at higher energies than implied 
by the fit form. The other is that, if one attempts the association 1/xmin Eqm, then 
eq. (|]) would have predicted 6 = 2e = 0.16, i.e. a much steeper slope. The alternative 
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association 1/xmin oc Eq^^ would give nicer agreement, but we iterate that it may be the 
less plausible ansatz. 

In summary, we have presented a very simpleminded model of the longitudinal mo- 
mentum, transverse coordinates and colour structure of the proton, in order to be able 
to estimate colour screening effects at small p±. Many of the details turn out to be less 
critical, but there are two key assumptions. The first is that the proton, by virtue of 
being a colour singlet, should couple to a soft gluon with reduced strength. The second 
is that, at higher energies, the proton can effectively be resolved into more partons that, 
in principle, can interact to give scatterings at a given p±_ scale. We show that these two 
assumptions lead to the conclusion that the perturbative ansatz for the cross section is 
only valid above some p^min scale that increases with energy. If this increase is fitted 
to a form p^mm oc Eq^, any value k in the range 0.05 to 0.16 could be obtained, which 
certainly overlaps with an experimentally acceptable behaviour. However, one lesson is 
that the ansatz in eq. may tend to overestimate the rate of increase at large energies, 
and thereby e.g. lead to an underestimate of multiplicities at LHC energies. Conversely, 
the study of minimum bias events at LHC may help improve the understanding of very 
low-p_L physics processes. 
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